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Abstract 

The connection of the 'time' evolution of the eigenstates of the refiectionless potentials 
of the Lax hierarchy to the more general case of the 'time' evolution of the eigenstates of the 
Schroedinger equation for potentials with non-vanishing reflection coefficients is explored. 
A new hierarchy of functions satisfying 'time' dependent equations is established. 



1 Introduction 

The Korteweg-deVries (KdV) equation was first discovered in the study of water waves. The 
KdV and the related equations with higher order nonlinearity, which are members of the Lax 
hierarchy (Lax 1968), have played a fundamental role in the study of nonlinear systems because 
they simulate many physical systems (Scott et al 1973), admit many conservation laws and the 
multi-soliton solutions can be given in analytic form. The KdV equation and its generalisation 
Kadomtsev-Petviashvli (KP) equation have also played an important role in pure mathematics 
because of their connection to algebraic curves, Jacobian varieties, vector bundles on curves, 
Schur polynomials and infinite dimensional Grassmannians (Mulase 1984). 

The connection between N-soliton solutions of the KdV equation 

dU d 3 U aTT du 

and refiectionless potentials with N boundstates in non-relativistic Quantum Mechanics is well 
known (Kay and Moses 1956, Gardner et al 1967, Scott et al 1973 ,Thacker et al 1978). The 
KdV equation is the m = 1 member of the Lax hierarchy of equations (Lax 1968, Caudrey et al 
1976) defined by 

dU 8L m 



dt m dx 



(2) 



where [Lf\ satisfy 



, p& a f)TT\ f>T 

L ° = U > (ft? -«> Tx - 2^) Lj . x = i = L2,.. )TO (3) 
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and t m is the 'time' parameter of the m th member of the hierarchy. The solutions to the non- 
linear equations of the Lax hierarchy may be used as potentials in linear Schroedinger equations 
and their spectral properties may be studied by solving (using units in which mass =1/2 and 
h=l) 

H = + U , H %l) k — E k ip k . (4) 

The N-soliton solution of eq. (2) may be viewed as a reflectionless potential U which supports 
N boundstates of the Hamiltonian operator H. Lax (1968) has shown that when the 'time' 
evolution of U is governed by eqs. (2) and (3), and the 'time' evolution of the eigenstates tp k is 
governed by 

dip 

B m ip k = ijr^- , ipk{x,t) = cxp(-iB m t) tp k (x,0) , (5) 
where the operator B m satisfies the commutator relation 

[Bm > H] = — ST = t aiZ = t aiZ' (6) 

then the eigenvalues E k of H are time independent and the eigenstates remain normalised, but 
the normalisation constants of the boundstates and the reflection and tranmission coefficients 
acquire a time dependence. For the case m = 1, which leads to the third order KdV, the explicit 
form of Bi is given by 

B ( 4 ^ 6U ^ 3^^^ (7) 

1 V dx 3 dx dx ) 

If the potential evolves in t\ according to eq. (1) and the eigenstates evolve in t\ according 
to eqs. (5) and (7) then the eigenvalues [E k = —7^] in eq. (4) are independent of t\, but the 
normalisation constants [C k ] which determine the behaviour of [ip k ] as x — >• ±00 and the reflection 
coefficient R for positive energies E — k 2 depend on t\ as given by (Scott et al 1973) 

C fc (7fc,ti) = C k ( lk ,0) exp(-4 7 |t 1 ), R(k,h) = R(k,0) exp (8 i k 3 h) . (8) 

Similar results hold for other members of the Lax hierarchy with appropriate energy dependent 
changes to C k (-y k ,t m ) and R(k,t m ). 

B m is a Hermitian operator which may be interpreted as the generator of 'time' evolution which 
propagates the potential U according to the KdV equation or another higher order non-linear 
equation arising from eqs. (2) and (3). This propagation of U is distinct from the usual time 
evolution of the Schroedinger eigenstates by the Hamiltonian H which propagates the parti- 
cle through a fixed potential. The hermiticity of B m ensures unitary 'time' evolution of the 
eigenstates of the Schroedinger equation. 

In this paper we develop methods by which is is possible to establish some general results which 
are valid in general and the Lax hierarchy is a special case of these general results. It will be 
shown that it is possible to construct a hierarchy of inter-connected functions which satisfy an 
evolution equation in parameter space which is different from eq. (2) but reduces eq. (2) for the 
case of the reflectionless potentials in a suitably taken limit. 

In section 2 of the paper we outline the steps involved in a known procedure for starting from a 
potential which supports no eigenstates but has non-vanishing reflection coefficient to constuct a 
potential with boundstates at specified energies and scattering properties related to those of the 
original potential in a well defined manner. In section 3 we study how a parametric evolution 



2 



different from eqs. (2) and (3) may be realised and establish some general equations. In section 4 
we discuss an example of a parametric evolution equation which may be related to the evolution 
of a physically meaningful quantity. 

When referring to 'time' in the rest of this paper we mean a parameter of the type which appears 
in evolution equations such as the KdV equation for the propagation of the potential. We do not 
refer to the time coordinate t which is canonically conjugate to the energy in Quantum Mechanics 
and appears in the evolution equation generated by the Hamiltonian. With this clarification we 
drop the quotation marks when referring to 'time' in the rest of the discussion. 



2 Construction of potentials with non-vanishing reflection 
coefficients 

The algorithm for the construction of a reflectionless potential U with TV boundstates starting 
from free particle states (Thackcr et al 1978) is a special case of a more general problem of 
starting from the solutions in a potential Uo which supports no boundstates and satisfies Uq — > 
as x — > ±oo and then finding a potential with boundstates at specific energies (Abraham and 
Moses 1980, Sukumar 1986, 1987, Baye 1987, Baye and Sparenberg 1994). Let <f>k be a solution 
of the Schroedinger equation for the potential Uq at energy E k = —7^ given by 

d 2 

Ha = +U , H <f> k = -ft 4>k (9) 

which goes to exponentially as x — > 00 but grows exponentially as z ^ -00 and therefore 
unnormalisable. A set of such solutions [cj)k] at energies [— 7I, k = 1,2, ..AT], may be used to 
define a matrix [A] with elements 

A u =6 kl - [ <p k {y) <pi{y) dy , k,l = 1,2, ...N . (10) 

The solutions [ipi] to the system of linear equations 

N 

Y,A M V; = 4>k , fc = l,2,..,A (11) 

2=1 

may be used to define 

N N N N N „ . „ 

W(x) = -5>(*)lb0«0 = -EE^" 1 ]^ = EE = oT m det A (12) 

1=1 1=1 fc=i 1=1 fe=i 

Eq. (11) may be written in the form 

N „ x 

Mx) = <f>k{x) + E / Mv)Mv) d V M*) , k = l,2,...,N (13) 
1=1 J °° 
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and can be differentiated twice and rearranged to give 

N „ r N 



N „t o N N 



;=i 1/00 i=i 1=1 

(14) 

Using the relation 

and eqs. (12) and (13) it can be established that 



(H - 2— + 7fe 2 )^(x) = J2J Mv)Mv)dy(H - 2— + 7?)^) 
which may be brought to the form 

( H o 2— + 7?)^ = 0. (17) 



(16) 



We can therefore infer that, if det A ^ 0, then 

ffjv = #0 + U N - U (18) 
dW d 2 

U N - U = - 2 — = -2— IndctA (19) 

(tfjv + 7?)^ = 0, 1= 1,2, ...,7V (20) 

(i.e.) [?/>z] are solutions of the Schroedinger equation in the potential Un for the energy — jf. 
Using eqs. (10) and (11) it can be shown that [ipi] satisfy boundstate boundary conditions at 
x — > ±00 establishing [ipi] as true boundstatcs of Un- 

It can also be shown that eq. (13) may be used to define tp(E,x) for positive energies E = k 2 
by the replacements tpk — > i>{E) and <pk — > <i>{E). With these replacements eq. (16) is still valid 
and hence, if eq. (20) holds, then 

il>(E,x) = <j>(E,x) + ]T / HE,y) My) dy M*) (21) 
1=1 J °° 

(H N -E)il>{E,x) = Y, / <K E >V) Mv) dy (H N + jf)M^) = (22) 
1=1 J °° 

(i.e) ip{E, x) is a solution of the Schroedinger equation for the potential Un at energy E. It may 
then be shown that for positive energies E — k 2 , the reflection and transmission coefficients of 
the two potentials U n and Uq are related by 

R N (k) = Ro(k) fl^rU, T n (k) = To(fc) {-)» f[^4l ■ (23) 
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A relation satisfied by the Wronskians of the two sets of solutions [fa] and [ipk] may be established 
by differentiating eq. (11), multiplying by fa and summing over k. These algebraic manipulations 
yield 

1=1 

The differentiation of eq. (12) gives 



9^ , %\ aw i 

(=i 

We can then establish that 

AT 



These results are valid for a general Uo and are the generalisations of known results for the case 
U Q = 0. 

This procedure for adding N boundstates to a potential Uo may be reversed so that starting 
from a potential with N boundstates in Un we can find a potential with no boundstates in the 
form 

d 2 

U = U N - 2 — IndctB (28) 

PX PX 

Bu = S M + / ipk(y) ipl(y) dy = I tp k (y) ipi(y) dy (29) 

J oo J — oo 

and the two sets of functions <j> and are now related by 

N N 

E B «<& = ^ > & = EC- 8-1 ]**^ ( 3 °) 
i=i 

so that 



fe=i 



JV N N N N 



i=l i=l fe=l 2=1 fe=l 



Comparison of the two sets of equations corresponding to the addition and the removal bound- 
states shows that the matrix [B] is the inverse of the matrix [A] which means that in terms of 
the kernel defined by 

N 

K{x,y) = EV-fc(z) fa(y) (32) 
k=i 
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eqs. (11) and (30) may also be given in the form 

N N 



(33) 



Mx) - E^Lfc^ = Y, Blk(t)k = ^ + / My)K(y,x)dy 
fc=i fc=i 

N N x 

4>i(x) = ^[-B" 1 ]^ = ^AfcV-fc = ^(s) - / K(x,y) Mv) dy . (34) 

I 1 1 1 '-'OO 



fc=l fe=l 
Eq. (19) may be expressed in terms of K(x,x) in the form 

dW d 

u n ~ U = -2 — = 2 g^K(x,x) . (35) 

3 Evolution in parameter space 

The potential Un constructed by adding boundstates to Uq may be viewed as the potential at 
time t = 0. Uq(x,0) and Un{x, 0) may now be used to construct time dependent potentials 
Uo(x,t) and J7jv(x,t) by letting both [7 and Un satisfy one of the non-linear time evolution 
equations of the Lax hierarchy, given by eqs. (2) and (3), with the same value of m, and 
letting the eigenstates and [ipk] evolve according to the appropriate set of eqs. (5) and (6). 
Under these circumstances the boundstatc eigenvalues remain unchanged while the boundstate 
normalisation constants and the reflection coefficients for positive energies for both potentials 
acquire a time dependence as discussed in section 1. However these changes will be the same for 
both potentials and the reflection coefficient and transmission coefficients for the two potentials 
will still be related in the manner indicated by eq. (23). Therefore, if we now eliminate the N 
boundstates of UN{x,t) by the procedure outlined in the last section the resulting potential will 
acquire a time dependence exactly of the form of Uo(x,t). Only for the case of reflectionless 
category of potentials, for which the reflection coefficient vanishes for all t, Uo(x,t) is zero for 
all t and U$ may be viewed as being time independent. This identification allows us to examine 
a different kind of parametric evolution in which Uq{x) is a potential which in general has non- 
vanishing reflection coefficient, but remains independent of time, and the energies [Eh] at which 
boundstates are added are also independent of time. 

Now suppose that a dependence of on a parameter t can be introduced through a set of 
functions [«&(£)] so that 

^ = a k (t) 4> k , <j>k(x,t) = cxp^ a k (i) dtj </> fc (x,0) (36) 

then <ftk(x,t) can still satisfy eq. (9) and Uo may be chosen to be independent of t. The 
corresponding t evolution of [ipk(x,t)} may be examined by differentiating eq. (11) with respect 
to t to get 

'dA M dtp, 



i=i 

Since 



— A M = (a fc + a t ) (A kl - 6 kl ) (38) 
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we can simplify eq. (37) using eq. (11) to the form 

N „ . 



1=1 

Hence the time evolution equation of tpi is governed by 

N 

- -I- n,lh,\ = 9 

v dt 



f^ti 4- a ,Vi) = 2 ]T [A- 1 ]^^ = 2 ]T B ;fe a fe Vfe (40) 



fc=i fe=i 

AT 



r x / \ 

= 2 Mv,t)^(Mv,t)Mt)Mx,t)) *v ■ ( 41 ) 
- 7 - 00 fe=i 

Multiplication of eq. (39) by tpk, summation over k and use of eqs. (11) and (36) then leads to 
the relation 

i=i fc=i 
which when combined with eq. (19) leads to 

f) r) r) N 

-(U N - U ) = j-U N =4- ^(s.t) . (43) 



• T/izts we /mue derived an evolution equation for U (x, t) which is valid for arbitrary choices of 
values of the functions [otk{t)] and depends on a sum over the boundstate probabilities with weight 
factors [afc]. 

Since [tpk] are eigenfunctions normalised to unity, integration of eq. (42) between — oo and oo 
leads to the relation 

d N 
£**->-oo — IndetA = 2Va t (44) 

when the result from eq. (10) that as x — > oo det A — > 1 is used. 

Using the Schroedinger equation satisfied by the eigenstates ipk it can be shown that the proba- 
bility density ipp. associated with the eigenstates of any potential satisfy 

(s-^- 1 f)*-«5«' l - 1 'W <«' 

The weighted sums 

N 

Qj = -4^/3 fe (2 7fe ) 2j ' $ , .7-0,1,2,.. , (46) 
k=l 

where [flk] can take arbitrary values, therefore, satisfy 

( 93 ATI 9 2 9Un \ O dQ * +1 (47) 

_ 4L/w Tx 2 ^x~) Qj ■ (47) 

If the time evolution functions in eq. (36) are chosen to be [ak] = {f3k{2"fk) 2 ^], then in terms of 
the members of the [Q] hierarchy, eq. (43) becomes 
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Thus eqs. (36), (10), (11) (18-20) and (46)-(48) generalise the equations for the Lax hierar- 
chy. 

For the case Uq = 0, the 7V-soliton potential Un = L where L is the m = member of the Lax 
hierarchy defined by 

N 

L m (x,t) = -2 ]T(2 7fe ) 2m+1 tf(x,t) , m = 0,1,2,... . (49) 
k=i 

Hence for the special case of the reflectionless potentials Un of the Lax hierarchy, constructed 
from Uo = 0, one of the weighted sums in eq. (46), (viz.) Qo, arising from the time independent 
choice of [/3fc(t)] = [7k], can be identified with the iV-soliton potential Un- All other [Qj] can be 
expressed in terms of Un and its derivatives using eq. (47) and eq. (43) becomes 

dU N dL m 

~dT = ~ ~dx~ (50) 

leading to the non-linear time evolution equation for the potential given by eq. (2). We have 
shown that for the case of Uo = the new parametric evolution discussed in this section becomes 
identical to the evolution equation of the Lax hierarchy given by eq. (2). This identity does not 
extend to potentials with non- vanishing reflection coefficients if we retain the choice of [(3k] = [jk] 
and in general a new non-linear evolution equation arises which differs in essential respects from 
eqs. (2) and (3). 

We have studied how far it is possible to go in formulating parametric evolution equations, 
for potentials with non-vanishing reflection coefficients, which differ from the time evolution 
equations of the Lax hierarchy . There exists the possibility that in the general case, for a 
suitable choice of the functions [«&(£)], which define the time development of the solutions in a 
time independent Uo, one of the weighted sums [Qj] can be related to a physically meaningful 
function. We examine such a possibility next. 



4 Green's function hierarchy and time evolution 

For the time independent choice 

«fc(t) = &(*) = 4(7 / + E) , E = - 7 2 , 7 [ 7fe ], (51) 

it is possible to relate Qo m eq. (46) to a Green's function by the following procedure. A 
Green's function G(x, y; t, E) satisfying suitable boundary conditions may be defined as a solution 
of 

(-^2 + U{x,t) - E)G{x,y;t,E) =S(x-y). (52) 

It is well known that G(x, y; t, E) may be expressed in terms of the complete set of eigenfunctions 
and eigenvalues of the Schroedinger equation for the potential U(x, t) (Morse and Feshbach 1953). 
For a general U(x,t) the complete set of eigenfunctions include eigenstates for both discrete and 
continuous eigenvalues and G may be represented by 

G(x, y;t,E) = - £ M^)My,t) J J&Lp^ x , t)m ^ t) dk (53) 
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where f(k) is the spectral density function for the potential U for positive energy E — k 2 . An 
alternative representation of G can be given in terms of £1 and £2, the two linearly independent 
solutions of the Schroedinger equation for U at energy E, 

{-^2 +U(x,t) - E)^ 2 { Xl t) = (54) 

chosen such that £i(x,t) — > as x — >• —00 and £2(^1 1)) — > as x — >• 00. The Wronskian of the 
two solutions at the same energy E defined by 

^^to-to^ (55) 

does not depend on x as a consequence of the differential equation satisfied by £1 and £2- Then 
G may be given in the form 

G(x,y;t,E) = - (56) 

where x< (x>) is the smaller (larger) of (x,y). By following the same method as that used to 
prove eq.(45) it may be established that 

(^-4(U-E)--2-)G(x,x;t,E)=0. (57) 
If the potential U is a confining potential with only a discrete spectrum then 

Qo = -E^4 = G(x,x;t,E) = (58 ) 
and G(x,x;t,E) satisfies the sumrule 

/oo ^ 
G(x,x;t,E)dx = -^-2-^ . (59) 

Eq. (43) now becomes 

-[/(x,t) = - —G(x,x;t,E) (60) 

Even if the potential is not a confining potential and has scattering states, the members of the 
hierarchy of functions defined by the discrete sums 

Gj(x,y;t,E) = - ^ - -^—^ , 3 = 0,1,2,.. (61) 

satisfy 

{^-< U »- E )Tx- 2 ^) Gj+l{x ^ E) = q^^x^E) • (62) 
Using the orthonormality of [ipk] it is also possible to establish the integral expression 

/oo 
Gjixw^E) G (x u y;t,E) dx x . (63) 
-OO 
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Using eq. (43) it may be seen that the time dependent potential constructed from eqs. (36), 
(10) and (19) for the choice [at — (47^ +4-E) -3-1 ] with eigenstates given by eq. (11) and Gj 
given by eq. (61) satisfies 

g-U N (x,t) = -—Gj(x,x;t,E) (64) 

• This result is valid in general and this time evolution equation for U in terms of the [G] hierarchy 
exists for all potentials including the special case of Uq = 0. In this new evolution equation Gj 
plays the same role as Lj does in the Lax hierarchy. It is possible that other choices of may 
exist for which the weighted sum ^ ctfcV'fc ma y represent a useful physical quantity. 

In this paper we have shown that the methods used to establish the time evolution of the 
reflcctionless potentials of the KdV equation can be extended to derive a parametric evolution 
equation for general potentials which is diffrent from that of the Lax-KdV hierarchy. We have 
identified a generalised Lax hierarchy of functions and derived a new type of implicitly non-linear 
time evolution equation for general potentials which produces a (unitary?) time evolution of the 
eigenstates of the Schroedinger equation for the potential. 
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